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Abstract 



O, We analyze the neutral spectrum of the holographic quantum Hall fluid de- 

^ scribed by the D2-D8' model. As expected for a quantum Hall state, we find 

I— —I the system to be stable and gapped and that typically the lowest excitation 

P^ mode is a magneto-roton. In addition, we find magneto-rotons in higher modes 

I> as well. We show that these magneto-rotons are direct consequences of level 

_^ crossings between vector and scalar modes. 
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1 Introduction 

The quantum Hall (QH) effect describes the behavior of electrons in a two-dimensional 
conductor in a transverse magnetic fieldjj When the filling fraction, defined as the 
ratio of charge density to the magnetic field in units of ^^ and denoted by z/, takes 
integer (integer QH effect) and certain rational (fractional QH effect) values, the 
electron fluid enters a QH phase with many striking properties. Most characteristic 
of these, of course, is the response to an applied electric fleld; the conductivity in the 
direction of the fleld vanishes, while the conductivity in the orthogonal direction, or 

2 

Hall conductivity, is exactly ^z^- 

The fluctuations of QH fluids exhibit many interesting properties which have been 
observed experimentally. In the QH state, a mass gap is dynamically generated for 
both the charged and neutral excitations. The charged fluctuations are quasiparticles, 
which in the fractional QH effect are fractionally charged and which obey anyon 
statistics. 

The lowest neutral excitations of quantum Hall fluids are magneto-rotons. These 
are collective excitations whose minimum energy is at zero group velocity but nonzero 
momentum. Rotons were flrst observed in superfluids, as local minima in the phonon 
dispersion at large momentum. The magneto-rotons in quantum Hall fluids occur 



^For a review of the quantum Hall effect, see u\. 



as global minima of the energy at nonzero momentum and have been detected in 
many experiments [2}]4]. There are several theoretical many-body and field theoretic 
treatments [5]-[7J. 

The important features of the QH effect are universal, in that they do not depend 
significantly on the details of the experiment, such as, for example, the crystal struc- 
ture of the material. This kind of generic phenomenon makes a promising candidate 
for holographic modeling. 

Our program is to investigate a class of top-down holographic models of strongly- 
coupled (2-|-l)-dimensional fermions exhibiting quantum Hall behaviorjj We have so 
far studied two similar brane constructions, the D3-D7' model of the fractional QH 



effect [9] and the D2-D8' model of the integer QH effect 10 . In both of these models, 
the QH state corresponds holographically to an embedding of the flavor brane which 
does not enter the black hole horizon. This feature of the embedding leads directly 
to a dynamically-induced mass gap and the expected conductivities for a QH fluid. 

The neutral fluctuations of the holographic QH state of the D3-D7' model were 
studied recently in 12 ,^ That investigation was limited to zero temperature. As 
expected, the spectrum was found to be gapped, indicating that the fermion fluid is 
incompressible. A tachyon was found in the branch of solutions previously identified 
as thermodynamically unstable fo] , but otherwise the QH state was stable. The lowest 
mode was shown to become rotonic in the vicinity of the crossing of the lowest scalar 
mode and the lowest vector mode. 

In this paper, we perform a similar fiuctuation analysis of the QH state of the 
D2-D8' system. We investigate the neutral bosonic spectrum and improve upon the 
D3-D7' study by considering nonzero temperature. As expected, we again find a 
mass gap, with a tachyon only in the thermodynamically unstable solution. We also 
show that, unlike in the D3-D7' model, the speed of sound is independent of mode 
number and weakly dependent on charge density. We find that, over a wide range 
of parameters, the lowest mode is a magneto-roton, and in addition, magneto-rotons 
appear in higher modes as well. We provide an analytical argument showing how 
these magneto-rotons are in fact direct consequences of level crossings between scalar 
and vector modes. 

In the next section. Section [2| we review the D2-D8' model, its action and so- 
lutions. We then set up the analysis of the fiuctuations and describe our numerical 
methods in Section [3j Section |4] contains our results, the spectrum and dispersion 
of the modes. We investigate the relationship between the observed level crossings 
and the rotonic dispersions in Section |5| Finally, we conclude in Section [6] with a 
discussion and open questions. 



^The QH effect has also been studied using several other types of holographic models [sl. 

■^An analogous fluctuation analysis of the ungapped, metallic state of the D3-D7' system was 



conducted in 
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2 Review of the D2-D8' model 

The D2-D8' system fTo] consists of a probe D8-brane in the near-horizon background 
of A^ thermal D2-branes, such that 4^ND = 6 in the flat-space limit. We work in 
the limit where A^ ^ 1 and QsN ^ 1. Supersymmetry is completely broken, and 
internal flux is needed for stability. At weak coupling, the low-energy spectrum of 
bifundamental strings in a H^ND = 6 system contains only charged fermions and no 
charged bosons; as a result, the holographic field theory dual is (2+l)-dimensional 
SYM coupled to A^ species of charged fermions. Similar constructions have been used, 
for example, in the D3-D7 \TT\ and D3-D7' [9] models. 

2.1 Set-up 

The ten-dimensional, near-horizon D2-brane metric is 

L~'^dslQ = u^ [-f{u)dt^ + dx^ + dy^) + WM — — + u^dVt^ 

where L^ = Gn'^ggNl^ and the thermal factor f{u) = 1 — (^) . There is a black hole 
horizon at m = ut, and the corresponding Hawking temperature is T = -^Urj. . The 
metric on the S^ is given by 

dVL% = di)"^ + sin^ i) [dOl + sin^ 9id(j)l) + cos^ ^ (c^C^ + sin^ ^dOl + sin^ i sin^ 612^02) > 

(2) 
where if) ranges from to 7r/2, ^, ^i, and 62 range from to vr, and (/)i and 02 range 
from to 27r. The background dilaton and RR 5-form potential are 

e* ^ «.(^)* (3) 

C(5) = c{ij)L^dVLs2 A rffis3 , (4) 

where c{;ip) = | (sin?/) — | sin(3^) — :j^ sin(50')). 

We will work in dimensionless (lowercase) coordinates, which are related to the 
physical dimensionful (uppercase) coordinates as: x^ = X'^/L, u = U/L. In general, 
lowercase Latin letters will denote dimensionless quantities, and uppercase letters will 
denote the corresponding physical quantities. 

The probe D8-brane fills the spacetime directions t, x, y, and the radial direction 
u and wraps an S*^ x S^ fibered over an interval in the internal S^. The embedding 
in the -^-direction is a function of u. The induced metric on the D8-brane is then 

L^'^dsj)^ = u^ {-fdf + dx'^ + dy'^) +u~^^ (- + u'^i)''^] du"^ 

+u~ ^ sin^ip dQl + u~^ cos^ ip dfll , (5) 



where ' = du- 

We would like to consider a boundary system with nonzero charge density in a 
background magnetic field. In the bulk, this is accomplished by turning on compo- 
nents of the worldvolume gauge field A^; specifically, a constant field strength in the 
spatial directions corresponds to a transverse background magnetic field, and nonzero 
charge density is dual to a radially varying electric field: 

27ra'F^y = h (6) 

2na'Fuo = aQ^u) . (7) 

In addition, we include a magnetic field wrapping the internal S"^ in order to stabilize 
the embedding: 

27ia'Fe^^^=L%sm(f)i . (8) 

Furthermore, we choose to work in radial gauge, where Au = 0. 
The DBI action for the D8-brane is then 



Sdbi = -fi8 I d\e-'^\/-det{g^, + 27Ta'Ff,,) (9) 

(10) 



-Af f duJ cos^ijJ{l + fu^ij'^ - a'^) {h^u + sin^^) M + 



h 



u^ 



where M = STr^TgfsL^ and v^ is the dimensionless volume of spacetime. As a result 
of the charge density and the magnetic field, the Chern-Simons action has a nonzero 
term: 

Scs = - — {27ia'y f C^^^ AFAF = Af f duc{ij)ha'o . (11) 

2.2 Background solution 



The D8-brane action (10) and (11) is independent of Oq, so the equation of motion 



for uq can be integrated once, giving 



where we have defined 



(r 2 
l + -]a', = d-hc{^) = d{u) , (12) 



9W = I = (13) 

:i + fu^i^'^ - a'^) (1 + S) 



and where d is the constant of integration and d{u) is the radial displacement field. 
In terms of the boundary theory, d is the total charge density, while d{u) is the charge 
due to sources located in the bulk at radial positions below u. 



We can solve for Oq in terms of d, obtaining 



1 + fu^i/j'^ 
d? + u^ cos6 i> [h^u + sin^ i)) (l + g) 



% = d\l ^^ . . . , /,n . ■ . ,^ /. . fe2^ • (14) 



We then rewrite (13) as 



/ /d2 + M5cos6^(l + g) (6% + sinV) 
l+^V 1 + fu^ij'' 






The ^/^ equation of motion, written in terms of d, is 

du I u^g I 1 H — ^ ) V^' ) ~ ^'^^^ V' sin ip (2 cos^ ^|J sin^ -^/^ — 3 (jj^u + sin"^ '^) ) 

f hd n , , 

-St^^^cos^^. (16) 



^1 + 



,5 



We can solve ( 16 ) at large u and find that t/^ — )■ in the UV. We identify the coefficient 
of the leading term in the asymptotic expansion of ip as the fermion mass and the 
coefficient of the subleading term as the chiral condensate: 

^ = - - 4 + ^(«"') • (17) 

There are two classes of D8-brane embeddings with different IR behavior: black 
hole (BH) embeddings, which enter the horizon at ut, and Minkowski (MN) embed- 
dings, which smoothly cap off at some uq > ut as the wrapped 5''^ shrinks to zero 
size. As argued in flOj, the MN embeddings holographically reproduce the properties 
of a quantum Hall fluid, while the BH solutions exhibit metallic behavior. In this 
paper, we will focus on fluctuations of the MN embedding, reserving the study of the 
quasinormal modes of the BH solutions for a future work. 

The QH features of the MN solution are direct consequences of the boundary 



conditions at the tip, iI){uq) = 7t/2. We found in 10 that regularity of the gauge 



flelds requires that the ratio of the charge density d and the magnetic fleld h be flxed: 

^=c(./2) = H. (18) 

The flUing fraction per fermion species u is deflned by the raticQof the physical charge 
density D = ^ '^'^J d to the physical magnetic fleld H = g^, divided by the number 
of species A^. We flnd that 



^The filling fraction v is given in units of the Dirac flux quantum ^^^, which is just ^^^ in units 
where e = h= 1. 



The MN embedding is therefore dual to an integer QH fluid. 



We found evidence in 10 to support this identification. By avoiding the horizon, 
the MN embedding generates a dynamically-induced mass gap for charged excitations, 
given by uq. In Section |4], we will find a mass gap for neutral fluctuations as well. 
In addition, as expected in a QH state, the Hall conductivity is given by the filling 
fraction, a^y = ^, while the longitudinal conductivity vanishes, a^x = 0. 

2.3 Isotropic coordinates 

To investigate the fluctuations of the D8-brane, it is much more convenient to work 
in Cartesian coordinates rather than the polar coordinates {u,ip). This is simple to 
accomplish at zero temperature; however, at nonzero temperature, because of the 
thermal factor /, we first need to convert to a new, isotropic radial coordinate r 
where 

\ 5/2 n TT 




-•''■^"Uj +1/1:^1 -1- (20) 

The horizon is located at r = 1, and the boundary is at r = cxd. For convenience, we 
define 

^ = i^ . ^ = i^. (21) 

This allows the thermal factor to be written as 

S-{g . (22) 

and the ten-dimensional metric takes the form 

L-Hsl, = {urrf'P (-tdf + dx' + dyA + ^y^^-^ {dr' + r'dnl) . (23) 

We can now introduce Cartesian coordinates on the (r, '?/')-plane: 

p = rsinip (24) 

R = r cosip , (25) 

where p and R give the sizes of the S^ and S^, respectively. We parametrize the 
D8-brane embedding as p{R), and the radial gauge field is 27ia'FjiQ = a'^i^R), where 
the prime now denotes partial derivatives with respect to R. The induced metric on 
the D8-brane is then 



L-^dsls = {uTrf^^^(-^df + dx^ + dyA 



1 



1/2^5/2^1/5 



^{{i + p'^)dR^ + p^dnl + R^dnl) . (26) 



The D8-brane DBI action in these coordinates is 



S 



DBI 



U / dR u^tR^0^/^, 



\ 



1 + P'2 - 






h^ 



y5^5^2 



P^ 



^^7-5^2/5 



+ 62 



and the Chern-Simons action is 



Scs = M I dR c{p)ha'Q , 



(27) 



(2J 



5p3 p5 



where c(p) - 3^ 

The background MN embedding caps off smoothly at r = ro as the S^ collapses, 
which implies the IR boundary condition p{R = 0) = tq. We still choose a radial 



gauge, which is now Aj^ = 0. Solving the equations of motion derived from (27) and 



(28), we find solutions for the embedding p{R) and radial gauge field ao(-R). The 
fermion mass is related to the asymptotic embedding by 



p(oo) = 2^/5 



r, m 



Ut 



(29) 



The mass as a function of tq for several values of d is shown in Fig. [TJ 




''0 



Figure 1: The fermion mass m plotted as a function of ro for ut = 0.1 and 6=1. 
The three curves correspond to different values of rf: rf = 0.1 solid blue, d = 1 dashed 
black, and d = 10 dotted red. Note that, embeddings that seemingly yield a negative 
mass are actually unstable due to self-intersections (Tol. 



3 Fluctuations 

3.1 Expansion 

We now allow for fluctuations of the following worldvolume fields: 

p = p(R) + e6pit,x,y,R) (30) 

do = ao{R) + e6aQ(t,x,y,R) (31) 

a^r = e6aa:{t,x,y,R) (32) 

tty = hx + e6ay(t,x,y,R) (33) 

ttR = e5aR{t,x,y,R) , (34) 



with e a small parameter and where a^ = 2^^m- ^^ ^'^ ^'^^ consider fluctuations 
in the internal spheres or dependence on the internal S"^ or S^ coordinates. The 
determinant det(G' + 2TTa'F) in the DBl action is now 



-{uTr)'/% e{^-6d,) e{^-6ay) a', + e{6a', - 6dn) 



-e 



13 ^ dx ^' ^ dy 

^-Sa,) {urrf/^P -h + ed^ySa,] e{6a',-^] 



(35) 



--a',-e{5a',-5dn) -e(K - ^) "^K " ^) ^7$^ 



where dots denote derivatives with respect to time and primes are derivatives with 
respect to R. The Chern-Simons action becomes 

Scs = A/" / dR c{p + e6p)e^ M ^ + 5a'o - 5a^ j I- + d[Jay] j 

The on-shell Lagrangian can then be expanded to second order in e: 

C = Co + eCi + e^C2 + ... . (37) 

Here Cq corresponds to the background Lagrangian, and £i vanishes by the back- 
ground equations of motion. We derive the equations of motion for 5p, 5ao, 5ax, and 
5ay from £2- However, since they are very long expressions, we will not reproduce 
them here. In addition, since we are working in the radial gauge or = 0, the equation 
of motion for da^i becomes a constraint to be imposed on the fluctuations. 
We make the following wavelike ansatz for the fluctuations: 

6pit,x,R) = (5p(i?)e-*"*+'^" (38) 

6a^{t,x,R) = (5a^(i?)e-*'^*+*'^% (39) 

9 



where, because of the rotational symmetry in the (x, ?/)-plane, we can choose the 
momentum to be in the x-direction without loss of generality. 

We must be careful to ensure we consider physical fluctuations rather than gauge 
artifacts. To ensure this, we work with the fluctuation of the longitudinal electric 
field: 

5ex = k6do + ujda^ (40) 

which is gauge invariant. As long as we start with initial conditions in radial gauge 
at some R which satisfy the constraint given by the equation of motion for 5aR, the 
constraint and gauge condition will be satisfied by the solution at all R. 

3.2 Numerical methodology 

Our goal is to solve the equations of motion for the fluctuations and find the normal- 
izable modes. The condition for normalizability is that the fluctuations vanish at the 
boundary: 

5ex{oo) = Sdy{oo) = 6p{oo) = . (41) 

Furthermore, the embedding of the D8-brane must be smooth at the tip, which implies 

6p\0) = , (42) 

and the gauge fields must also be nonsingular, giving 

6e'M = K(0) = . (43) 

In practice, to numerically solve the coupled equations for the fluctuations at 
fixed values of u and k, we vary the initial conditions at the tip and shoot out to the 
boundary to look for solutions which are normalizable. We must solve three coupled, 
linear, second-order differential equations for three independent functions. 



An efficient technique for finding such solutions is the determinant method 14 , 15 
Rather than attempting to scan, for each u and k, the full three-dimensional space 
of initial conditions, we choose a basis of linearly-independent initial conditions, 

i6p{0),6eM, S^ym = {(1, 0, 0), (0, 1, 1), (0, 1,-1)} (44) 

for which we numerically solve the equations of motion. We then compute the fol- 
lowing determinant of boundary values: 



5p^ 5p^^ 5p^^^ 
det I 6ei 6ei^ 5ei^^ 



(45) 



where the superscripts refer to the basis elements of (44). If the determinant vanishes 



there is at least one linear combination of the basis vectors for which 6p{oo), Sex{oo), 

10 



and Say{oo) simultaneously vanish, signaling the existence of a normalizable mode. 
We repeat this procedure for various values of u and k; for fc = 0, a scan in u yields 
the spectrum, and allowing k > enables us to compute the dispersions. 

For the figures presented in Section |4| we have chosen particular values of the 
parameters ut, rn, d, h, and b for illustrative purposes. As we will be studying 
exclusively MN embeddings, we will work with the parameter d and always set h = ^d. 
The scaling u — )■ Am, ut -^ ^ut, m — )■ Am, d — )■ A^^'^d, h — )■ A^/^/i, and b — )■ A~^/^6 
is a symmetry of the equations of motion, and so there is, in fact, only a three- 
dimensional parameter space of MN solutions. Many of the qualitative results we 
found appear to be independent of these choices. However, as an exhaustive numerical 
search is essentially impossible, we can not conclusively rule out the possibility of 
different behavior in particular corners of the parameter space. 

^0 




Figure 2: The mass ojq of the eight lowest modes plotted against ut, for d = 10, b = 1 
and m = 0.1. Blue solid curves are vector modes and red dashed curves are scalars. 
There are three mode crossings: modes one and two cross at Ut ~ 1.25, modes four 
and five cross at ut ~ 1.50, and modes seven and eight cross at ut ~ 1.54. The 
maximum temperature for the MN embedding is ut ~ 1.548. 



4 Numerical results 

4.1 Spectrum 

We compute the spectrum by finding the normalizable solutions at fc = 0. The 
masses uq = uj{k = 0) of the lowest eight modes as a function of temperature for 
representative values of the other parameters are shown in Fig. |2j As expected for a 
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Figure 3: The mass gap as a function of d for ut = 0.1, b = 1, and m = 0.1. The 
blue circles show cjq for the lowest vector mode, and the red squares for lowest scalar 
mode. The two levels cross at rf ^ 70. The magenta diamonds, which lie slightly 
below the Uq of the lower mode, show the minimum energy u^, in the region d > 2.36 
where the lowest mode is a magneto-roton. The mass gap is shown by the curve, 
which is solid blue for d < 2.36 and dashed magenta for d > 2.36. 



QH fluid, we see a mass gap; the size of the gap, shown Fig. [3} is given by the energy 
of the lowest mode. 

The mass squared is positive for all modes, indicating perturbative stability, with 



one exception. In 10 we found in the vicinity of the high-temperature phase tran- 



sition that there are two branches of MN solutions, and the MN embedding with 
smaller tq was found to be thermodynamically unstable. The fluctuation analysis 
shows that this small-ro MN solution is perturbatively unstable, as well. 

Fig. E^ shows Uq for the lowest two modes for both large-rg and small-ro branches 
as functions of temperaturejj The lowest mode of the unstable branch (dashed red 
curve) is tachyonic. The two branches meet at ut ~ 0.67, just as the lowest mode of 
the stable solution (solid red curve) becomes massless. 

Another feature we find at fc = is a decoupling, also observed in 12 , of the scalar 



Sp fluctuation and the vector {6ex, Sdy) fluctuations. At low temperature, the lowest 
mode is a vector, followed by a scalar. Subsequently, every third mode is another 
scalar. At nonzero k, the equations of motion are fully coupled, but we can continue 
to employ their k = labels. 

From Fig. [2| it is evident that all of the masses decrease with increasing temper- 
ature, though at different rates. In particular, the masses of scalar modes decrease 
faster than the masses of vector modes, leading to periodic levels crossings. As the 
temperature is raised, the crossing of the two lowest modes, shown in Figs. [2] and |4| 
occurs first. At higher temperatures, the fourth and fifth modes cross and then the 
seventh and eighth modes, as shown in Fig. [2j While there is the potential for a cross- 



^For ease of comparison, we have chosen the same parameters for Fig. |4] as for Figs. 7 and 
in Tol. 
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Figure 4: The mass squared uj^ of the two lowest modes plotted against mt, for 
d = h = 1 and m = 0.1. The solid curves are the fluctuations of the stable MN 
solution; blue is a vector and red is a scalar. The dashed lines are the fluctuations of 
the unstable solution. The dashed red curve has uJq < 0, with a positive imaginary 
part, and is therefore tachyonic. 

ing at every third mode, for a given choice of parameters, a finite number of crossings 
will occur before the maximum temperature is reached and the perturbatively stable 
branch of MN embeddings terminates. 

We have described the spectrum as a function of increasing temperature at certain 
fixed values of the charge density (and magnetic field). However, the behavior is 
qualitatively the same if we hold the temperature fixed and decrease the density. 

4.2 Dispersion 

We now consider the dispersion of these modes for nozero k. In most cases, the modes 
have massive dispersion relations, well-approximated by 



u 



^0 + c^k 



21,2 



(46) 



where Cg is the speed of sound. However, near the crossing of a scalar and vector 
mode, the lower of the two modes develops a rotonic dispersion, with a minimum 
u^: < Uq located at A;* > 0, which is well approximated by the form 



uj= \/ujI + cl{k - kj' 



(47) 



as is shown in Fig. [5j However, near k = 0, the dispersion instead takes the form 
cj ~ Wo — 0{k'^). This deviation is a result of the level crossing and will be addressed 
in Section m 
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Figure 5: The dispersion of the lowest mode aX d = 1, ut = 0.61, 6 = 1, and m = 0.1. 
The black dots are the numerical result, the solid red curve is the fit to the rotonic 

1.360 and fc^, 
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form (47) with u^, 



fit to (54) with ujq 



1.0867, and the dashed blue curve is the small- /c 
0.0377. 



0.5 



0.4 



0.3 



0.2 



0.1 



0.1 0.2 0.5 1.0 2.0 5.0 10.020.0 



d 



Figure 6: The speed of sound c^ for the lowest mode plotted against the charge density 
(i in a log scale for ut = 0.1, 6 = 1, and m = 0.1. Note that for small d, the mode is 
massive, while for rf > 2.2, it is rotonic. 



We compute the speed of sound by fitting the large k numerical data to (46) and 



(47). We find that c^, defined by the large k asymptotics of the dispersion, is fully 



independent of the mode number and whether the mode is massive or rotonic. The 
speed of sound is weakly dependent on the charge density; it increases slightly as d 
is reduced, as shown in Fig. [6j 



14 



We will first focus on the magneto-roton appearing in the lowest mode; the 
magneto-rotons in higher modes are qualitatively similar. We show in Fig. [7] the 
region in the {d, ut) parameter space where this magneto-roton exists. For small 
values of the charge, d < 2.4, the magneto-roton appears for a narrow range of 
temperatures, while for large charge, the magneto-roton exists all the way to zero 
temperature. Although the upper bound on the magneto-roton range is quite close 
to the upper bound on the existence of the MN embedding, there is a small gap 
between the two. 




Figure 7: The region in the (c?, Mj')-plane in which the lowest mode is a magneto- 
roton, for 6 = 1 and m = 0.1, bounded by the blue solid curves. The black dashed 
curve is the location of the first-order phase transition, and the dotted black curve is 
where the MN embedding solution ceases to exist. Although the region near d = Q 
is numerically difficult to investigate, we have verified that the magneto-roton exists 
at least down to d = 0.1, where the upper and lower bounds are at most 5ut ~ 0.02 
apart. 



Specializing for definiteness to (i = 1, Fig.|8]shows the dispersions of the lowest two 
modes for four different temperatures, and Fig. [9] plots the location and energy of the 
magneto-roton minimum, k* and w*, as functions of the temperature. This magneto- 
roton exists in a finite range, 0.52 ^ uj- ^ 0.665 around the level crossing, which for 
these parameters occurs at ut ~ 0.62. Note that there is a small but nonzero distance 
between the upper bound and the maximum temperature, ut ~ 0.67. Right at the 
level crossing, the magneto-roton minimum reaches a maximum k^,. 

The magneto-rotons observed in the fourth and higher modes are qualitatively 
much the same as the magneto-roton in the lowest mode. They occur near level 
crossings between scalar and vector modes, although the regions of parameter space 
over which they exist are much smaller. Fig. 10 is the analog of Fig. [8] for the fourth 
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Figure 8: The dispersions uj vs. k of the lowest two modes, the scalar in red and 
the vector in blue, are shown at four temperatures, ut = 0.51,0.60,0.62, and 0.64 
with d = b = 1 and m = 0.1. The upper left panel shows ut = 0.51; the modes 
both have the standard massive dispersion. In the upper right, where ut = 0.60, the 
vector mode has become rotonic. At ut = 0.62, in the lower left, the two modes are 
degenerate at fc = 0, so the scalar and vector modes can not be distinguished; we 
denote this by plotting both in black. In the lower right, at ut = 0.64, the scalar is 
now a magneto-roton. 
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Figure 9: The location k^ of the magneto-roton minimum (left) and the minimum of 
the magneto-roton energy u^ (right) plotted as functions of temperature Ut- Both 
figures are ioi d = b = 1 and m = 0.1. 
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Figure 10: The dispersions u vs. k of the lowest five modes, scalars in blue and vectors 
in red, are shown at four temperatures, ut = 1-3, 1.45, 1.50, and 1.52 with d = 10, 
b = 1 and m = 0.1. In the upper left panel, at ut = 1.3 there is a magneto-roton 
only in the lowest mode. In the upper right, where ut = 1.45, the fourth mode, a 
vector, has become rotonic as well. At ut = 1.50, in the lower left, the fourth and 
fifth modes are degenerate at A; = and so are shown as black, with the fourth mode 
rotonic. In the lower right, where Ut = 1.52, the scalar is now the fourth mode and 
is a magneto-roton. 



mode, showing the rotonic dispersion at several temperatures. As is evident from 
Fig. [2] and Fig. 10, the crossings of higher levels, and thus the higher magneto-rotons, 
occur at higher temperatures. As the temperature is increased, we generally find a 
narrow region where two magneto-rotons coexist. 



5 Level crossing and magneto-rotons 

Let us perform a brief analytical analysis to clarify the relationship between the level 
crossings and consequent rotonic dispersion relations. The solutions to the fiuctuation 
equations of motion define a linear function which maps the IR boundary conditions 
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to the UV ones. We may define its 3 x 3 matrix representation as 

5hy{00) \=M\ 5hy{{]) . (48) 

5p(oo) / V 5p(0) / 

The matrix depends on the external parameters ut and d as well as on uj and k. 

As pointed out above, for /c = the off-diagonal elements of M mixing the scalar 
5p and vector {5ex,5ay) fluctuations vanish. At fc = 0, we can diagonalize M by 
choosing an appropriate basis for the vector fluctuations, 

5e± = a 6ex ± i(3 Shy , (49) 

where a and /3 are real coefficients which must be determined numerically. Notice 
that this change of basis is not a unitary transformation. 

Consider further a choice of ut and d such that a scalar mode and a vector mode 
have almost equal energies, say cj-t = t^o ± 5, where 5/ujq ^ 1, and the plus sign 
corresponds to the vector mode. Let us fix the coefficients a and (5 such that M is 
diagonal dXuj = uj+ and A; = 0, with 5e- corresponding to the normalizable fiuctuation 
mode. We order the new basis vectors as (5e+, 5e_,(5p). Then only Mn and M33 are 
nonzero in the new basis when evaluated at the point defined abovejj Expanding to 
the first order in 5 and k around this point gives 

(aii(a;-cj+) + cii 012(0; -a;+) \ / 613 \ 

-ai2{uj-u+) a22{uJ-u+) + ° 623 U • 

033(0; -o;_) / V ^31 &32 / 

(50) 
Notice that the linear fc-dependence only appears in the mixing between the scalar 
and vector modes. This refiects a discrete symmetry of the fiuctuation equations 



of motion. Indeed we can check that the Lagrangian densities of (35) and (36) are 
invariant up to O(e^) (but not at higher orders) under the transformation 

x^ -X ] y\-^-y; da^^y H- -6ax,y (51) 

which leads to the symmetry 

k H- —k ; 5dy H- —5dy ; 5ex ^-)■ —5ex (52) 

of the equations of motion for the fiuctuation. Also the matrix M must be covariant 
under this transformation, which explains its structure. 
For normalizable modes we require, as above, 

det M = . (53) 

^We could also let a and (3 depend on u such that at A; = 0, M would be diagonal for all w. 
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Calculating the energies of these modes from the approximate matrix, to first order 
in 6 and k, gives 



wi 



W2 



W3 



CUq- 



W0 + 



U+ 



6^ + ^^^F 



0-22033 



52 + ^^^k^ 

0220-33 
C11O22 



O11O22 + 0'12 



(54) 

(55) 
(56) 



The dispersions (54) and (55), with the values of aij and bij extracted from numerics. 



give the small-fc behavior of modes near a level crossing and match the behavior 
observed in Figs. |5| |8| and[l0j In particular, the lower mode, Ui, initially decreases 
with k. On the other hand, at large k, the energy will instead be increasing like c^/c; 
there must necessarily be a local minimum at some nonzero k. This magneto- roton 
minimum is therefore a direct consequence of the level crossing. 



Another feature of (54) and (55) which we find in Figs, ^and 10 is the change in 



fc-dependence at the point where the modes are degenerate. When 6 > 0, the first 
/c-dependent correction to the energy is quadratic: a; ~ w-t + O(fc^). However, when 
the modes are degenerate and S = 0, the dependence on the momentum is linear 
u ^ uq + 0{k) . The dependence on the expansion coefficients appears only trough 
the combination 623^32/(022033), which is invariant under arbitrary scalings of the 
basis vectors (5e_ and Sp. In order for the energies to remain real as k increases, this 
combination must be positive. While it is difficult to exclude the apparently spurious 
case of negative 623^32/(022033) by analytical methods, it did not show up in any of 
the level crossings that we studied numerically. Notice also that the result depends 
only on the mixing coefficients 632 and 623 between the two nearby modes. In fact, 
ignoring the other mode and simply diagonalizing the lower-right 2x2 subblock of 



M gives exactly (54) and (55). 



The energy of the additional vector mode 003 does not depend on k to linear order, 
as is expected for modes in the absence of level crossing. Since this mode does not lie 
within ~ (5 of the modes which are about to cross, the linear expansion gives at best 
a rough approximation to its energy. However, we may improve the precision of 003 
by adding higher order corrections in u in the calculation, which does not change the 
fc-dependence of M as it is protected by the discrete k — )■ —k symmetry. Therefore, 
linear terms in k will continue to be absent in the expression for u^. 

6 Discussion 

We continued in this paper our investigation of the D2-D8' model by analyzing the 
fluctuations of the QH state at nonzero temperature. We found that the neutral 
spectrum is stable and gapped, providing evidence that the MN embedding describes 
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an incompressible QH fluid. We also showed that the thermodynamically unstable 
phase is perturb at ively unstable, as well. 

As observed in real QH fluids, the lowest mode was found to have a rotonic 
dispersion, at least for much of the parameter space. In addition, we uncovered 
magneto-rotons in higher modes as well. We showed that all of these magneto-rotons 
are direct consequences of near degeneracies at zero momentum between a scalar 
mode and a vector mode. Experiments have detected two such closely spaced, long- 



wavelength modes 16 . In 17 , a theoretical treatment has been given in terms of a 



hydrodynamical model. A composite fermion model has also been proposed 18,19 
and identifles the second mode as a two-magneto-roton bound state. 

Many of our results for the D2-D8' model are qualitatively similar to those found 
in [12] for the D3-D7' system. While both models found magneto-rotons in the lowest 
mode, in the D3-D7' model a magneto-roton was only found in a narrow region in 
parameter space and no magneto-rotons were seen in higher modes. However, D3- 
D7' investigation was restricted to zero temperature, and it potentially could be that 
magneto-rotons are more common at higher temperatures. 

Of course, these types of i^ND = 6 brane constructions are far from perfect 
holographic models of QH systems, which is not surprising given that they contain 
in their construction many unphysical features, such as large A^. Furthermore, there 
is a parametrically large hierarchy between the mass gap for charged excitations and 
for neutral excitations, Uq ""^^^ j^neutrai ^ _\2/5^ which is absent in real quantum Hall 
systems. However, the quantum Hall effect is a very robust phenomenon in (2+1)- 
dimensional fermion systems, so we are hopeful that there is a large universality class 
of models which exhibit its characteristic properties. The close qualitative similarity 
between the D2-D8' and D3-D7' systems, we believe, suggests that these types of 
brane models lie within such a universality class. 

We restricted our attention in this paper to the normal modes of the MN em- 
beddings of the D8-brane; in a future work, we will address the properties of the 
quasi-normal modes of the BH embeddings, which holographically model a Fermi-like 
liquid. The analogous quasi-normal mode analysis of the D3-D7' model was per- 



formed in 13 . The neutral spectrum was found to be ungapped and to feature a 
zero-sound mode. Interestingly, at sufficiently high charge density, a tachyon appears 
at nonzero momentum, indicating an instability to a spatially-modulated, spin and 
charge density wave state. As both brane models contain similar Chern-Simons terms 
of the type which generate this sort of instability, we anticipate that it should appear 
in the D2-D8' system, as well. 
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